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Abstract 

Background: The development of ultrasound for use in dental tissues is hampered by the complex, multilayered 
nature of the teeth. The purpose of this preliminary study was to obtain the phase and group velocities associated 
with several directions of ultrasonic wave propagation in relation to the tooth structure, which would then lead to 
the determination of the elastic constants in dental hard tissue. Knowledge of these elastic constants can be used 
to feed back into numerical models (such as finite element) in order to simulate/predict ultrasonic wave 
propagation and behavior in the teeth. This will help to optimize ultrasonic protocols as potential noninvasive 
therapeutic tools for novel dental regenerative therapies. 

Methods: An extracted human second molar was used to determine time-of-f light information from A-scan 
signatures obtained at various angles of inclination and rotation using a scanning acoustic microscope at 10 MHz. 
Phase and group velocities and associated slowness curves were calculated in order to determine the independent 
elastic constants in the human teeth. 

Results: Results show that as the tooth was inclined at three azimuthal angles (0 jn = 0°, 15°, and 30°) and rotated 
from (Din = 0° to 360° in order to cover the whole perimeter of the tooth, slowness curves constructed from the 
computed phase and group velocities versus angle of rotation confirm the inhomogeneous and anisotropic nature 
of the tooth as indicated by the nonuniform appearance of uneven circular shape patterns of the measurements 
when compared to those produced in a control isotropic fused quartz sample. 

Conclusions: This study demonstrates that phase and group velocities of ultrasound as determined by acoustic 
microscopy change and are dependent on the direction of the tooth structure. Thus, these results confirm that the 
tooth is indeed a multilayered anisotropic structure underscoring that there is no single elastic constant sufficient to 
represent the complex structure of the tooth. The findings underline the importance to take into account these crucial 
characteristics in order to develop and optimize therapeutic as well as diagnostic applications of ultrasound in dental 
tissue repair, and further studies are warranted to analyze ultrasound transmission at various frequencies and intensities 
in different teeth to develop appropriate models for ultrasound biophysical behavior in dental tissues. 
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Background 

Ultrasonic waves such as those used routinely during 
medical examinations and prognoses propagate in a very 
complex manner. This level of complexity holds true 
even when these sound waves propagate in ideal, homo- 
geneous, and uniform media. Despite the fact that 
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diagnostic ultrasound techniques are based on well- 
established phenomena, there remains the area of 
dentistry where ultrasound has not played yet a very im- 
portant role. This is due to the fact that ultrasonic 
propagation in inhomogeneous/anisotropic media results 
in further degradations in acoustic field. Scattering ef- 
fects at the microscale level as well as at the multilayer 
interfaces are responsible for these degradations. 

In earlier diagnostic ultrasound experiments performed 
on the human teeth, an ultrasonic pulse-echo system was 
able to detect the enamel-dentin junction and the dentin- 
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pulp interface [1]. Ultrasonic imaging techniques have also 
been applied to other areas of dentistry, specifically pros- 
thodontics. A B-mode ultrasonic system for determining 
the thickness of the masticatory mucosa for denture con- 
struction has also been developed [2], and the use of the 
scanning acoustic microscope to image the elastic proper- 
ties of carious human dental enamel has been described 
[3]. Acoustic microscopy has been successful for imaging 
lesions in tooth enamel and detecting the enamel-dentin 
interface [4]. Finite element models [5,6] as well as PSpice 
(OrCAD™, Cadence Design Systems, Inc., San Jose, CA, 
USA) transmission line (T-line) models [7] used to under- 
stand ultrasonic wave propagation in the teeth have been 
established. In order to confirm predictions made by these 
latter finite element and T-line models, experimental data 
was gathered [8] for the purpose of assessing the viability 
of ultrasound to detect cavities and fractures in the 
extracted human teeth. The method showed great promise 
as a potential improvement upon current dental imaging 
systems. 

For years, the main applications of therapeutic ultra- 
sound in dentistry have been mostly limited to oral sur- 
face cleaning (i.e., removal of plaque and calculus) or 
root canal treatment [9]. In recent years, however, thera- 
peutic ultrasound studies performed on dental tissue 
have taken a different path [10,11]. The goal of this re- 
search is to develop novel ultrasound therapies that can 
be easily applied to the teeth to promote tissue healing, 
more specifically tooth regeneration. Preliminary results 
showed that this idea of repairing the teeth with ultra- 
sound has potential. The tooth is a hard mineralized 
structure with a living soft tissue core capable to re- 
spond to outside stimuli including physical forces. Re- 
search has demonstrated that ultrasound in the low 
frequency range, generally used in dentistry for dental 
scaling, stimulated gene expression and the production 
of tissue regulatory growth factors that are considered to 
be crucial for odontoblast differentiation, activity, and 
dentin repair [12,13]. More recently, it was demon- 
strated that low-frequency ultrasound stimulated differ- 
entiation and in vitro mineralization in odontoblast-like 
cultures [11]. These findings are novel and exciting and 
support the working hypothesis that ultrasound may be 
an effective noninvasive tool for novel dental regenera- 
tive therapies [14]. Interestingly, these findings were re- 
cently corroborated in in vitro tooth slice organ cultures 
where low-intensity pulsed ultrasound was shown to in- 
crease pre-dentin thickness and sub-odontoblast cell 
numbers [10,15]. 

However, to develop ultrasound as a therapeutic appli- 
cation, a thorough understanding is required regarding 
ultrasonic wave propagation and intensity distribution in 
teeth structures. The tooth is a highly anisotropic 
medium. This quality stems from two different criteria — 



layer consistency and structural component. The enamel 
portion of a tooth mainly consists of calcium hydroxy- 
apatite [16]; the dentin represents a collagenous mineral- 
ized substance much like the bone and contains tubular 
canals, and the tooth's central cavity (pulp) contains a 
loose connective tissue with various cell types, blood 
vessels, and nerves. Odontoblasts surround the dental 
pulp with cellular processes extending into the dentinal 
tubules which may also contain neurite endings. In 
addition, the root is enclosed by a thin layer of bone-like 
material called cementum, which is surrounded by a 
periodontal ligament. The latter tissue is highly cellular 
and contains bundles of thick collagenous fibers respon- 
sible for attaching the tooth to the alveolar jaw bone. So, 
as it may be appreciated, the tooth is composed of a very 
complex mesh of intricate layers, each having different 
chemical and structural characteristics than the one adja- 
cent to it. As a result of all these effects, ultrasonic waves 
propagating within a tooth would undergo a fair amount 
of mode conversion. This, in a way, is similar to the as- 
sessment of the structural integrity of anisotropic compo- 
nents used in aircrafts and submarines whereby elastic 
constants of the host material need to be determined. 

In this paper, we propose and test a proof-of-concept 
methodology based on previous studies [17-21] that would 
eventually lead to the determination of the elastic con- 
stants in dental hard tissue. The study used an extracted 
human tooth to determine phase and group velocities 
from measured time-of-flight (TOF) of an ultrasonic wave 
as it propagates along arbitrary directions through the 
tooth crown of the specimen using a pulse-echo setup at 
different inclination angles of incidence. 

Methods 

All studies were approved by and in compliance with the 
guidelines set by the institutional review board commit- 
tee at Hofstra University (Hempstead, NY, USA). 

Anisotropy consideration 

Assuming a tooth can be classified as a simple ortho- 
rhombic material, nine independent elastic constants 
need to be determined. These are represented by the 
matrix in Equation 1: 
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The elastic constants along the diagonal are defined in 
the <1 0 0>, <0 1 0>, and <0 0 1> directions. In other 
words, the waves that are traveling along these directions 
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are considered to be pure mode, and their velocities 
(i.e., longitudinal (V L ), shear horizontal (V S h)> and shear 
vertical (V sv )) are determined along these directions. Note 
that all diagonal elastic constants can be determined by 
propagation waves in symmetric directions. The waves 
traveling along the lower-symmetry directions are not 
pure modes. The behavior of elastic waves propagating in 
anisotropic materials [17] is governed by the following: 



C 13 = ±{[2pV 2 ~(C n + C 33 + 2C 5 5)] 2 -i(C33-C 11 ) 2 }-C 5 5, 



(?) 



C 12 = ±{[2pV 2 -Uc ll + C 33 + 2C 66 )] 2 --(C 22 -C ll ) 2 f-C 66 . 

1 4 



(8) 
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where C^i is the elastic constant, u t is the displacement, 
and p is the mass density of the material. It is assumed 
that a plane wave solution u b with wave vector /<, fre- 
quency (o, and amplitude A b of the form u t - A t e l(<kx ~ ^ 
exists and can be substituted back into the equation of 
motion (Equation 2) in order to obtain Kelvin- Chris toff el 
elastic stiffness equation as follows: 



(r ik -p V 2 8 ik )u k = 0, 



(3) 



where T ik - Ci^aiaj are Christoffel elastic stiffnesses 
[17], Fis the phase velocity, and a t and a y are the direc- 
tion cosines of propagating wave. A general solution to 
(Equation 3) gives relationships between the phase velo- 
cities Fand Cij for the wave propagation along arbitrary 
directions. As an example, in the <0 1 1> direction, we 
should have the following: 



C 55 + C 66 0 0 
0 C22 + C44 C23 + C44 

0 C23 + C44 C33 + C44 



(4) 



then, using the solution of Kelvin-Christoffel's equation, 
we get the following: 

pV 2 =]-(C 55 + C 66 ) 



or 



pV 2 =^{C 2 2 + C 33 ± [(C 22 -C 33 ) 2 + 4 (C 23 + C 44 ) 2 ]'} + C44. 

(5) 

If Cu is known and pV 2 is measured, then C 2 3 can be 
calculated as follows: 



C 23 = ±{[2pV 2 -^(C 22 + C 23 + 2C 44 )] 2 -i(C 33 -C 22 ) 2 }^-C 44 . 



(6) 



Similarly, using the velocities along the <1 0 1> and 
the <1 1 0> directions, one can obtain C 13 and C 12 as 
follows: 



The above description says that one can measure the 
wave velocities along the <1 0 0>, <0 1 0>, and <0 0 1> 
directions for the diagonal elements C n , C 22 > C 33 , C 44 , 
C 55 , and C 66 and measure the wave velocities along the 
<0 1 1>, <1 0 1>, and <1 1 0> directions for the off- 
diagonal elements C 23 , C 13 , and C 12 . One should appre- 
ciate that this is not a trivial procedure. In order to 
measure the velocities along the aforementioned direc- 
tions, four similar (if not exact) teeth specimens need to 
be prepared as follows: 

1. All surfaces are perpendicular to the symmetry 
direction in such a way that the length in each 
direction is at least five longitudinal wavelengths for 
proper measurement, and all surfaces are flat. 

2. Largest surface perpendicular to the <0 1 1> 
direction and upper and lower surfaces are flat. 

3. Largest surface perpendicular to the <1 0 1> 
direction and upper and lower surfaces are flat. 

4. Largest surface perpendicular to the <1 1 0> 
direction and upper and lower surfaces are flat. 

However, as mentioned earlier, this study is only a 
proof-of-concept and was performed on only one tooth 
sample. Future studies will consider the aforementioned 
implementation to measure the elastic constants from 
obtained phase velocities. 

Experimental setup 

It should be emphasized that teeth come in a variety of 
shapes and sizes (different morphologies— e.g., incisors 
and molars). Also, the teeth are so brittle that extreme 
care must be taken during the cutting process so that 
they do not get fragmented. In this study, an extracted 
adult human (age, sex, or race not known) second molar 
was tested using a phosphate buffered saline (PBS) solu- 
tion as the coupling medium. The tooth was sterilized 
with an autoclave system and was completely dry with 
no vital root. In fact, the latter was not needed in this 
study since the focus was on the determination of mater- 
ial properties associated with specifically the dentin layer 
as described below. Since the study deals with the deter- 
mination of physical tooth characteristics (wave veloci- 
ties and eventually elastic constants), both the top 
surface of the crown and the bottom cusp surface 
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surrounding the pulp chamber were flattened by grind- 
ing to maintain good wave propagation that was not 
hindered by uneven surfaces. This resulted in a minia- 
turized tooth size of d = 6.25 mm in length (as opposed 
to the 15-mm original length) and 10 mm in width 
(measured across the breadth of the crown area) as 
shown in Figure 1. However, some of the indentations 
due to natural occlusions (OC in Figure 1) were left. 
Also, one should keep in mind that due to the flattening 
process, the majority of the enamel layer has been lost, 
exposing only the dentin portion of the tooth. This in 
fact is advantageous since the major anisotropy/inhomo- 
geneity is present in the dentin due to the presence of 
collagenous mineralized substance and tubular canals, as 
mentioned earlier. 

The A-scan signals were acquired using a scan- 
ning acoustic microscope consisting of an ultrasonic 
pulser/receiver (Model 5800PR, Olympus-Panametrics, 
Waltham, MA, USA), an immersion transducer (Model 
V312, 10 MHz, 1.27 cm (0.5 in.) focal length, 0.635 cm 
(0.25 in.) aperture diameter, Olympus-Panametrics, 
Waltham, MA, USA), a FlexSCAN-C ultrasonic C-scan 
system (Sonix Inc., Springfield, VA, USA) with a tank 
containing PBS solution as the couplant, and a computer 
software system with a front-wall-follower gating cap- 
ability. A Plexiglas™ tooth holder (Figure 2) was ma- 
chined to hold the tooth under test in an upright 
position with the top of the crown facing the radiated 
acoustic field. A special rotating and inclining base was 
constructed in order for the ultrasonic beams to propa- 
gate at various angles in the tooth. Since the transducer 
could not rotate around the tooth or its axis, the device 
was made to rotate 360° around a pivot in the tank and 
enabled the determination of the TOF values at different 
angles (rotational 0 in and azimuthal/inclination @ in ). 
During data acquisition, it is important to note that the 
tooth was scanned with the transducer perpendicular to 
the tooth surface of the crown as shown in Figure 3. 




Figure 1 Flattened adult human second molar specimen used 

for study showing the natural occlusions (OC). 

\ J 



The transducer is designed for negative spike excita- 
tion. The maximum spike excitation voltage is 300 V 
with a fast rise time and short duration. Although nega- 
tive spike excitation is recommended, continuous wave 
or tone burst excitations may be used. It is helpful to 
note that the number of cycles per tone burst excitation 
depends on the duty cycle, which itself depends on sev- 
eral other factors such as impedance, total power, volt- 
age, and phase of the transducer. The spatial length of 
the pulse for the transducer used in this experiment has 
been determined to be about 0.2 us, which is small 
compared to the echoes obtained from the tooth layers 
(enamel, dentin, tubular canals, and pulp). The pulser/ 
receiver was set to send an ultrasonic pulse with an 
energy setting between 50 and 100 uj and a gain be- 
tween 20 and 40 dB to insure proper propagation in the 
tooth sample. The tooth (mounted in the fixture) was 
placed in the water tank in the focal plane within the 
transducer's cone of insonification. The FlexSCAN- 
C® system was used to examine the front-wall reflection 
off the crown and all subsurface layers so that the TOF 
can be adequately measured for speed determination. It 
is crucial to point out that the resulting A-scan is com- 
posed of multiple modes (longitudinal, edge longitudinal, 
edge shear, head, and surface) of the ultrasound wave 
transmitted as a result of mode conversion [8] . As an ex- 
ample, Figure 4 shows a typical RF A-mode scan (ampli- 
tude vs. time) from a normal unflattened tooth sample 




Figure 2 Plexiglas tooth holder showing the inclining O jn 
platform {PL) that holds the tooth (7). 9 in is controlled by a 
turning knob {KB) that controls a lever {LV) connected to {PL). A base 
pivot {PV) controls the angle of rotation (P in . 

V J 
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Figure 3 Experimental setup in the water tank. It shows the relative position of the ultrasound transducer (XD) and the tooth sample 
(77-/; shown in red circle) in the water tank (WT). 



similar to the one we are investigating. The boundary of 
each of the three major layers (enamel, dentin, and pulp) is 
clearly indicated in this figure. However, in our case, since 
the enamel layer has been removed due to the filing/flat- 
tening of the tooth, only the peaks corresponding to the 
dentin and the pulp were used for calculating the TOF. 

As noted earlier, the propagation of ultrasonic energy 
impinging on the anisotropic sample at arbitrary 
(orthogonal or oblique) angles results in three mode- 
converted waves. These can be observed by measuring 
the TOF along the energy path. The phase velocities 
needed for the determination of the elastic constants are 
a function of the propagation time. For an arbitrary 



angle of incidence (0 in and @ in ), the phase velocity [17] 
is given by the following: 

\fAt{G- m ,Q- m )\ 2 At(G in , 0 in )cos0 in / 1 \ 2 1 4 

VmM = [{ — d — ) ~ 2 — v^d — + \y~ w ) ' 

(9) 

where At is the one-way TOF, d is the thickness of the 
tooth (as shown in Figure 1), and V w is the phase vel- 
ocity in the PBS solution. Furthermore, using Snell's law, 
the angle of refraction 0 r can be determined from the 
following: 



< 
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Figure 4 A typical received ultrasound A-scan signal from a normal tooth and corresponding layers. (£, enamel; D, dentin; P, pulp) It 
shows how the TOF is determined. 
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cos@ r 
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v u 



V p (@ in , (Pin), 



(10) 



From Equations 9 and 10, the magnitude and direc- 
tions of phase velocities for all modes of propagation in- 
side the tooth can be determined for any arbitrary angle 
of incidence, @ in . Once this is done, an inversion routine 
based on the Newton-Rhapson method [17] can be 
implemented on the above characteristic polynomial to 
determine the elastic constants. Knowledge of these elas- 
tic constants can be used to feed back into models (such 
as finite element) in order to simulate/predict ultrasonic 
wave propagation and behavior in the teeth. 

Furthermore, once the angle of refraction @ r has been 
calculated for each phase velocity, one can determine 
the various group velocities V g which correspond to 
these respective phase velocities at the calculated angles 
@ r . These group velocities represent the actual family of 
longitudinal stress-wave velocities in those particular di- 
rections within the tooth. Group velocities can be deter- 
mined by projecting their vectors, on the surface of the 
slowness curves (as will be shown later), normal to the 
corresponding phase velocities. Using simple trigonom- 
etry, the group velocities can then be computed from 
the following: 



V, 



Vn 



COS0 r 



(11) 



Results and discussion 

Several tests were performed for obtaining TOF infor- 
mation using three azimuthal angles (angles of inclin- 
ation) set at @ in = 0°, 15°, and 30°. For each of these 
three angles, the platform was rotated from 0 in = 0° to 
about 360° in order to cover the whole perimeter of the 
tooth. The phase velocity in the PBS solution (V w ) is as- 
sumed to be similar to that in water (1,500 m/s). For 
each set of incident angle (0 in , @ in ), TOF and the phase 
velocity V p were determined. Polar plots of TOF and 
slowness curves from the computed phase velocities are 
then constructed. In order to confirm the experimental 
measurement system, TOF measurements as a function 
of 0 in and @ in were made on an isotropic, homogeneous 
square piece of fused quartz plate of 4-mm thickness. 
The polar plot of TOF values At ((& iw 0° -> 2n) f @ in ) 
determined from around the quartz plate is perfectly cir- 
cular for propagation of longitudinal wave, with a con- 
stant longitudinal velocity of about 5,714 m/s. This 
observation confirms the homogenous characteristic of 
the control quartz plate. Figure 5 shows the results for 
both TOF and longitudinal velocities obtained for the 
quartz sample at an angle of inclination (@ in = 0°). 
Figure 6 shows the TOF measurements for the tooth for 
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195 
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180 
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Figure 5 Panoramic polar plots of (a) one-way TOF (us) and (b) 
measured velocities (m/s) for quartz sample. The plots were 
generated at an angle of inclination 0 in = 0°. 



all three angles of inclination (@ in = 0°, 15°, and 30°), 
and Figure 7 illustrates the corresponding phase veloci- 
ties in the tooth for the same angles. 

The direction cosines of the phase vectors in the tooth 
can be determined from Equation 10. These are shown 
in Figure 8. 

The observations highlight polar/slowness curves 
(Figures 6 and 7) which confirm the inhomogeneous and 
anisotropic nature of the tooth. This is indicated by the 
nonuniform appearance of uneven circular shape patterns 
of both TOF and V p measurements when compared with 
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180 

Figure 6 Panoramic polar plot of the one way TOF (us) in tooth 
for all three angles of inclination, (a) 0 in = 0°, (b) 0 in = 15°, and 

(c) e in = 30°. 




180 

Figure 7 Corresponding panoramic polar plots of phase 
velocities (m/s) in tooth for all three angles of inclination. 

(a) 9 in = 0°, (b) 9 in = 15°, and (c) 9 in = 30°. 
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the perfect circular measurements made in the control 
isotropic fused quartz sample shown in Figure 5. This out- 
come has obvious implications for the development of 
ultrasound as a diagnostic tool, and in this respect, further 
analyses need to be conducted to determine the phase 
velocity of ultrasound at different frequencies. 

The results as shown in Figures 7 and 8 will help in the 
complete determination of the magnitude and directions 
of phase velocities for all three modes of wave propagation 
(quasi-longitudinal and two quasi-shear) inside the tooth 
for any arbitrary angle of incidence. These will then be 
used to determine the independent elastic constants of 
this orthotropic medium. Also, using Equation 11 in con- 
junction with the TOF and V p measurements shown in 
Figures 7 and 8, the group velocities can be calculated. 
Figure 9 shows the corresponding polar plots. 

As mentioned earlier, since the flattening process of 
the tooth caused the majority of the enamel layer to be 
lost, exposing only the dentin portion, the average group 
velocities attained were approximately 3,178, 3,088, and 
3,080 m/s for each of the three azimuthal angles (@ in = 
0°, 15°, and 30°), respectively. These values confirm what 
was reported in earlier studies [1,5-8]. 

Moreover, regarding the potential use of ultrasound for 
dental tissue repair, it is important to relate the physical 
behavior of ultrasound through the different layers of the 
teeth to those in a viable pulp cavity and relate these char- 
acteristics to localized ultrasound intensity in order to as- 
sess the biological effects of ultrasound in the pulp region. 
In other words, knowledge of the elastic constants can be 
used to feed back into numerical models (such as finite 
element) in order to simulate/predict ultrasonic wave 
propagation and behavior in the teeth. This will help 
optimize ultrasonic protocols as potential noninvasive 
tools for novel dental regenerative therapies. 

Conclusions 

This proof-of-concept study demonstrated that phase 
velocities (and corresponding group velocities) of ultra- 
sound as determined by acoustic microscopy change and 
are dependent on the direction of the tooth structure. 
The data underscore the anisotropic structure and 
demonstrate that no single elastic constant will be suffi- 
cient to represent the complex nature of the tooth. The 
findings have obvious implications for the development 
of diagnostic and therapeutic ultrasonic devices for 
dental tissues, and further correlation studies that will 
examine other samples are warranted to analyze ultra- 
sound transmission at various frequencies and intensities 
in different teeth to develop appropriate models for 
ultrasound biophysical behavior in dental tissues. 
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